Brane Decay of a (4+n)-Dimensional Rotating 
Black Hole. II: spin-1 particles 



M. Casals 1 , P. Kanti 2 and E. Winstanley 3 

1 School of Mathematical Sciences, University College Dublin, 
Belfield, Dublin 4, Ireland 

2 Department of Mathematical Sciences, University of Durham, 
Science Site, South Road, Durham DH1 3LE, United Kingdom 

3 Department of Applied Mathematics, The University of Sheffield, 
Hicks Building, Hounsfield Road, Sheffield S3 7RH, United Kingdom 

Abstract 

The present works complements and expands a previous one, focused on the emission of 
scalar fields by a (4 + n) -dimensional rotating black hole on the brane, by studying the 
emission of gauge fields on the brane from a similar black hole. A comprehensive analysis 
of the particle, energy and angular momentum emission rates is undertaken, for arbitrary 
angular momentum of the black hole and dimensionality of spacetime. Our analysis 
reveals the existence of a number of distinct features associated with the emission of spin- 
1 fields from a rotating black hole on the brane, such as the behaviour and magnitude of 
the different emission rates, the angular distribution of particles and energy, the relative 
enhancement compared to the scalar fields, and the magnitude of the superradiance effect. 
Apart from their theoretical interest, these features can comprise clear signatures of the 
emission of Hawking radiation from a brane-world black hole during its spin-down phase 
upon successful detection of this effect during an experiment. 



1 Introduction 



One of the most exciting consequences of the scenario of Large Extra Dimensions [1 
is that it opened the way for probing physics beyond the scale of quantum gravity by 
postulating the existence of a low fundamental gravitational scale M*. If this scale is 
sufficiently low, trans-planckian particle collisions could easily be achieved in various 
environments. The products of such collisions cannot be ordinary particles any more but 
rather heavy or extended objects (p-branes, string balls, string states, etc) arising in the 
context of the more fundamental theory valid above M*. As this ultimate fundamental 
theory should also include gravity, strong gravitational phenomena could also be observed 
at the quantum level, and the creation of tiny black holes as the result of trans-planckian 
collisions t 2j is highly likely. 

According to the theory with Large Extra Dimensions, all Standard Model fields 
are restricted to live on a (3+l)-dimensional brane, while gravity is allowed to propagate 
in the (4 + n)-dimensional bulk (for some early works, see The theory predicts the 
existence of n additional spacelike compact dimensions, all having - in the simplest case - 
the same size L. The collision of highly energetic particles on the brane may then produce 
a tiny black hole that is attached to our brane but, being a gravitational object, extends 
off our brane as well. In order to be able to treat the black hole as a classical object, it will 
be assumed that the black hole mass Mbh is considerably larger than the fundamental 
scale of gravity M*. Under the additional assumption that the horizon radius of the 
produced black hole is much smaller than the size of the extra compact dimensions L, this 
object can be considered to live in a (4 + n)-dimensional non-compact, empty 1 space; its 
properties differ from a 4-dimensional black hole with the same mass and strongly depend 
on the number n of additional dimensions present 

Small higher-dimensional black holes may be created during trans-planckian colli- 
sions either at ground-based colliders [S] or in high energy cosmic-ray interactions in the 
atmosphere of the earth [U] (for an extensive discussion of the phenomenological impli- 
cations and a more complete list of references, see the reviews |H1 El)- After a short 
'balding' phase, during which the black hole will shed all additional quantum numbers 
apart from its mass, angular momentum and charge, the more familiar Kerr-like phase 
commences, during which the black hole will lose mainly its angular momentum. Af- 
ter that, a Schwarzschild phase follows with the, now spherically-symmetric, black hole 
gradually losing its actual mass. The energy loss of the black hole takes place through 
the emission of Hawking radiation ^U] (and through superradiance during its Kerr-like 
phase) that is, at the same time, the most prominent signature of black hole creation. 
This consists of the emission of elementary particles both in the bulk and on the brane, 
and it will be characterized by a very distinct thermal spectrum. 

Due to the simplicity of the corresponding gravitational background that describes 
it, the Schwarzschild phase of the life of a small higher-dimensional black hole underwent 
an intensive and thorough study during recent years. The emission of Hawking radiation 
from such a black hole has been studied both analytically [TT] 1X21 113] and numerically |14j . 
These works have provided us with both analytical formulae and exact numerical results 

lr The brane self-energy can be naturally assumed to be of the order of the fundamental Planck scale M* 
in order to avoid a hierarchy problem. But since M* <C M^Hi its effect on the gravitational background 
is negligible. 
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for the emission rates, and demonstrated how the latter strongly depend on the existence 
of additional spacelike dimensions in nature. As was shown in detail in jT3], the number 
of particles and amount of energy emitted by the black hole per unit time was strongly 
enhanced by the presence of extra dimensions in the theory. Remarkably, it was found 
that the type of particles emitted also depends on how many dimensions exist in the bulk: 
scalar particles are predominantly emitted by a black hole living in a spacetime with a 
small number of (or no) extra dimensions, while gauge bosons are the preferred particles 
emitted by a black hole living in a spacetime with n > 5. The spectrum of Hawking 
radiation may be a source of valuable information on other fundamental parameters of 
the theory, such as the value of the cosmological constant [T3] and the string coupling 
associated with the presence of higher-derivative gravitational (Gauss-Bonnet) terms [16J. 
Other studies have focused on the effect of the mass of the emitted particles on the 
spectrum [H| HHj as well as the charge of the black hole (THj • 

Although initially almost totally ignored, the progress in studying the emission of 
Hawking radiation during the Kerr-like phase in the life of a small higher- dimensional 
black hole has been very rapid during the last year. Two early works [2HI E] focused on 
the special case of a 5-dimensional rotating black hole: the first one provided analytical 
formulae for the emission rates, and the second derived analytical results for the energy 
emission rate in the limit of low energy and low angular momentum of the black hole. 
The derivation of exact results valid for arbitrary values of the energy of the emitted 
particle and for arbitrary values of the angular momentum of the black hole required the 
numerical integration of both radial and angular parts of the equation of motion of the 
specific particle. The same task demands also the determination of the angular eigenvalue, 
that links radial and angular equations together, and which does not exist in closed form 
for arbitrary energy and angular momentum. The first exact numerical results for the 
emission of Hawking radiation from a (4 + n) -dimensional rotating black hole, in the form 
of scalar fields, as well as for the energy amplification due to the superradiance effect, 
were presented in [221- There it was shown that the energy emission is greatly enhanced 
by both the dimensionality of spacetime and the angular momentum of the black hole. 
In [2B1, results on the power (energy) spectrum, as well as its angular distribution in the 
5-dimensional case were presented, but once more the analysis relied on the assumption 
of low energy and low angular momentum. Results on the superradiance effect in the case 
of a higher- dimensional black hole were also presented in [21] and [22] , while the emission 
in the bulk was addressed in j2Hl I2Z] • 

Recently, a comprehensive study of the emission of Hawking radiation from a higher- 
dimensional black hole in its spin-down phase appeared in the literature [2E|- This work 
focused on the emission of scalar fields on the brane, and presented exact results for the 
fluxes of particles, energy and angular momentum, for arbitrary values of the dimen- 
sionality of spacetime, angular momentum of the black hole, and energy of the emitted 
particle. Also, the angular distribution of the emitted radiation - a distinctive signature 
of emission from a rotating black hole - was investigated: it was demonstrated that the 
uniform distribution of particles in the case of a spherically-symmetric black hole is now 
replaced by a clearly oriented one that is transverse to the rotation axis of the black hole. 

As the techniques used in the aforementioned work were strongly spin- dependent, 
the analysis for the emission of higher-spin fields was left for subsequent works. Here, we 
undertake the effort to fulfil this task for the case of gauge fields. As these fields are by 
definition restricted to live on the brane, only the brane emission channel will be studied. 
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To this end, we will use semiclassical techniques that were originally developed for the 
case of 4- dimensional black holes and are here generalized to cover the case of higher- 
dimensional black hole line-elements projected onto a brane. As these techniques were 
developed particularly for the case of spin-1 particles, we are forced to leave the study of 
fermions for another work. 

The outline of our paper is as follows: in Section 2, we present the theoretical 
framework for our analysis that includes the form of the gravitational background, and 
the derivation of the radial and angular part of the spin-1 equation of motion on the 
brane. In Section 3, we review, and generalize where necessary, the principles leading to 
the expressions for the transmission (or absorption) coefficient and of the various Hawking 
radiation emission fluxes. The techniques used to numerically integrate both the angular 
and radial parts of the spin-1 equation of motion - to derive the angular eigenvalue, the 
spin-weighted spheroidal harmonics, and the radial function - are described in detail in 
Section 4. The readers who are not interested in numerical techniques can move directly 
to Section 5, where our results are presented. Subsection 5.1 presents our results for the 
transmission coefficient and its behaviour in terms of the number of extra dimensions 
and angular momentum of the black hole. Subsections 5.2, 5.3 and 5.4 are devoted 
to the computation of the particle, energy and angular momentum fluxes, respectively: 
their dependence on both the dimensionality of spacetime and angular momentum of the 
black hole is investigated and, in the case of the first two types of fluxes, their angular 
distribution is also studied. Subsection 5.5 presents the total emissivities from a higher- 
dimensional rotating black hole on the brane in terms of the angular momentum, number 
of additional dimensions and azimuthal angle. The final subsection 5.6 addresses the 
issue of superradiance, and the energy amplification is computed in terms of the angular 
momentum and dimensionality of spacetime. We finish with a summary of our results 
and conclusions, in Section 6. 

2 Gauge Field in a Brane-Induced Rotating Black 
Hole Background 

As is well known, the generalization of the 4-dimensional, Kerr black hole background 
to a rotating, uncharged black hole living in a (4 + n)-dimensional spacetime is given by 
the Myers-Perry solution j2H] ■ Under the assumption that all ordinary matter - including 
ourselves - is restricted to live on a brane embedded in this background, we are interested 
in the type of line-element that describes the geometry on our brane. This follows by 
projecting the (4 + n)-dimensional Myers-Perry background onto the brane, and is found 
to have the form [2j 





(1) 



where 



A = r 2 + a 2 



•n— 1 



£ = r 2 + a 2 cos 2 9 . 



(2) 



r 



The mass and angular momentum (transverse to the r<^-plane) of the black hole are then 
given by 

(w + 2)A n+2 2 

M B h = // , J = — — M BH a , (3) 

167rG n + 2 

with G being the (4 + n) -dimensional Newton's constant, and A n+2 the area of an (n + 2)- 
dimensional unit sphere, given by 

27r (n+3)/2 

An+2 = r[(n + 3)/2] • (4) 

Above, we have assumed that the (4 + n)-dimensional black hole was created by the 
collision of two highly energetic particles moving on our 3-brane. Assuming that the 
thickness of our brane is much smaller than the size L of the extra dimensions, the two 
particles will have a non-zero impact parameter only along our brane; the produced black 
hole will therefore acquire only one non-zero angular parameter about an axis (the one 
transverse to the rcp-plsnae) in our brane. 

The black hole horizon is given by solving the equation A(r) = 0, which, for n > 1, 
leads to a unique solution given by 

where we have defined a* = a/r^. We should note here that while for n = and n = 1, 
there is a maximum possible value of a that guarantees the existence of a real solution to 
the equation A = 0, for n > 1 there is no fundamental upper bound on a and a horizon 
always exists. An upper bound can nevertheless be imposed on the angular momentum 
parameter of the black hole by demanding the creation itself of the black hole from the 
collision of the two particles. The maximum value of the impact parameter between the 
two particles that can lead to the creation of a black hole was found to be 



1 + 



n + 2 



i 



yU^+D , (6) 



an analytic expression that is in very good agreement with the numerical results produced 
in [2 c. Then, by writing J = bMsn 1 2 [21], for the angular momentum of the black hole, 
and using Eq. (jHJ) and the second of Eq. (|3*|). we obtain 

«r x = ^ • (7) 



We now proceed to derive the equation of motion of a gauge field in the induced-on- 
the-brane gravitational background The quantization of an electromagnetic field on 
a rotating black hole background is considerably more involved than for a non-rotating 
black hole. While there is an extensive literature on the quantization of a scalar field on a 
Kerr black hole spacetime, rather less work has been done on spin-1 fields. For example, 
the first comprehensive study of the stress-energy tensor on the entire geometry has only 
been recently undertaken jHU]. However, the formalism itself has been sufficiently well 
developed for our purposes for Kerr black holes (SHI EH E2 ESI EH EE1 EH EH EH] and, 
here, we shall extend this to the more general background metric In this section, we 
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will only outline the key features of the theory, referring the interested reader to [30J (and 
references therein) for further details. 

We use the Newman- Penrose formalism jSHj to derive the equations for the electro- 
magnetic field components, written in terms of the Newman- Penrose scalars: 

<?i-i = f/ m ; 0o = {Fi n + Fm*m) /2; 4>+i = Fm*n, (8) 

where {/, n, m, m*} are the null basis vectors. Here, we follow the notation of [HUHHIIEIj 
and use the subscript h on <ft, taking the values (0, ±1), to describe the helicity of the 
spin-1 field - note that both helicities (h = ±1) will contribute to the Hawking radiation. 
We should emphasize that, in our terminology, helicity is in effect a label for the field 
modes, and we do not strictly link this to any particular physical field configuration. This 
latter point is explored in more detail in [3UJII2I- 

The tetrad we use is the Kinnersley tetrad [321; which readily generalizes to the 
metric ((TJ with A given by Eq. (J2j) rather than the usual Kerr metric form. 2 The analysis 
of Teukolsky |3TJ E2] showed that the partial differential equations for 0_ x and p~ 2 + i 
[where p — — l/(r — ia cos 8) is one of the spin coefficients] are in this case separable. It is 
straightforward to show that this result also extends to our more general situation. The 
Maxwell scalars (jHJ) can be decomposed into a Fourier mode sum 

/+oo +oo +1 

dUJ 2^ h a lmu)P h4>lmu > (9) 

■°° e=\h\ m=-e p=±i 

where h a emujp are the coefficients of the Fourier series. The parameter i labels the eigen- 
values of the angular Teukolsky equation (see Eq. (|T^|) below). The parameter P is 
summed in (jUJ) over the values +1 and —1, corresponding to the two linearly independent 
polarization states of the electromagnetic potential. An electromagnetic potential mode 
of polarization P is an eigenstate of the parity operator 

V: (9,(f>) -> (n-9,0 + n). (10) 

with eigenvalue P. The decomposition of the potential into eigenstates of the parity op- 
erator V is the natural choice due to the invariance of the metric (0) under this operation. 
The Fourier modes in Q are then separated according to 



±Wtmw — ±1 



N emj p (1±1) Tl J? to (r) T iZ irmi {6, V ) e-™ , (11) 



where ±iN^ muJ is a suitable normalization constant, ^iRtmuj{ r ) is a radial function and 
^iZ(.muj{9, f) is the angular function. The equation for the Maxwell scalar O is not 
separable; instead, this scalar can be written in terms of the radial and angular functions 
for h — ±1. The answer is given explicitly in [30*1 133| 135] f° r Kerr black holes (n = 0), but 
here we do not need it further. The Fourier modes of the electromagnetic potential can 
be reconstructed from the Maxwell scalars ±i<fiirnu> - explicit formulae for this are given 
in |HnilH3]) again for Kerr black holes. 

The angular function hZimuiQ^) is defined as 

h Z emLU (6, <p) = 1 j= h S emw (6) , (12) 
v 2ir 



2 See [7||2] for details of the Newman-Penrose coefficients in this modified tetrad. 



where h St m u 1 {9) is the spin- weighted spheroidal harmonic EHl H5] l4T| Bgj 139], 
Introducing the new variable x = cos 6, the latter satisfies the equation 



jL ( (i _ ^ + fl2w v - 2^ - (m + f 

ax \ ax/ 1 — x 2 



hSe m cu(x) = 0. (13) 



In the above, hAi muJ is the separation constant, usually written as 



hMmw = h^imu - a 2 u 2 + 2mau, (14) 

with h^imuj regarded as the eigenvalue of the equation (fT3]). The spin- weighted spheroidal 
harmonics satisfy the following symmetry relations: 

-hSe muJ (^ — 6) 

-hSt-m-uiO), (15) 

and are normalized according to 





("1 


l+rn 


Sl mU }{9) = 


("1 


) h 




("1 


h+m 



d9mx9 h Si mu (9) = l. (16) 

The solution of the angular equation (|T5|l is discussed in Section E| 

The radial function h,Rimu>{ r )i on t ne other hand, is found to satisfy the following 



equation 

dr \ dr 
For h = +1, the potential hV = + \V is given by |2T] 



/^h±(^^^L \ +b y h n linj: = 0. (17) 



if 2 _ iXA' 

+i^ = ^ + 4ia;r+(A"-2)- +1 Aw, (18) 

where we have defined: K = (r 2 + a 2 ) ui — am. Throughout this work, we use a prime (') 
to denote a derivative with respect to the radial coordinate r. For both Kerr (n = 0) and 
5-dimensional Kerr-like black holes of the form given in Eq. (0), the term A" — 2 vanishes 
identically. The h — +1 radial equation (JT7jl with the potential (JTSJl fits the pattern of the 
master equation found for various spins in ■ However, there is a rather surprising 
result for helicity h = —1, when the radial function -l-ff&n^r) still satisfies (|17j) but the 
potential _iV has a different form: 

K 2 + iKA' . . 

_iV = — 4zu;r - „iA w . (19) 

Note that the A" — 2 term is absent in this case. Thus, while the master equation in [HE] 
describes the helicity h = +1 radial functions, it does not describe the helicity h — — 1 
radial functions. We suspect that something similar may happen with spin- 1/2 fields as 
well. The solution of the radial equation (JT7J) is also described in detail in Section HJ 



fi 



3 Hawking radiation in the form of gauge fields 



For an evaporating black hole on the brane, the relevant quantum state is the ("past") 
Unruh vacuum \U~). There are some subtleties in the construction of this state for an 
electromagnetic field on a rotating black hole background which are described in detail 
in |201 - T ne field is expanded in terms of the usual "in" and "up" modes, for which the 
radial functions hRimuij) have the asymptotic forms [30J: 



hRfrmb) ~ \ :„ -,„„ „-iun ■,„ „„ f „+iwr» . . (20) 



E)in,tra \-h a - 








pin,inc e~ iu3T * 
h^trruo r 


D in,ref e +ibJr * f„ 
I h n tmuj r l+2h V 


■* +oo) 




pup.inc +iwr* 


j_ pup,ref a -h „-i&r* 


(r - 


-> r h ) 


pup.tra e +i " r * 
hft-tmw r l+2h 




(r - 


-> +oo 



fc-EfoL ~ S „+io,r« . . ( 21 ) 



In the above, we have used the abbreviation 

Cj = uj — mQ, (22) 

with Q the horizon angular velocity 

fl= f1 A , (23) 
(1 + a*)r ft 

and the usual tortoise coordinate r* defined through the relation 

dr* r 2 + a 2 



dr A(r) 



(24) 



The coefficients h R^™, fc#*£?> ^Z'. appearing in Eqs. (EDD 

and (|21j) are arbitrary integration constants. The correct expression for the unrenormal- 
ized expectation value of the stress-energy tensor in the state \U~) is given in terms of a 
mode sum as follows [3U] : 

(U-MU-) = i£(T^cothg) 

<mP V>/U \ +/ 

+ ^ {T, v [ h <pf muj , h <PZ«] + (-1)^ {T,u {hVL^ h <t>tL] ) }) , (25) 
where the quantity [h^moj, h^lmJ * s the classical electromagnetic stress-energy tensor 

Tp,v [hfamu, h^tmJ = {-l<f>tmu -\^\mu n ^ n v + 2 </W> O^L ['(A* n «0 + m (/^ m ^)] 

++l<f>£mu +\4>\ m J'plv — ^0<Plmu) -l4>lmu J n {lJ, rri u) 

+complex conjugate, (26) 



given in terms of the tetrad basis vectors. In (|25|). the variable $ has been defined so 
that (—1)^ is equal to —1 if one index of the component of the stress tensor is 8 and the 
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other one is not, and it is equal to +1 otherwise. The key feature of the unrenormalized 
(U~\T^ U \U~) (J25)) is the term involving the parity operator V (|1U|). This ensures that the 
stress-energy tensor as a whole is invariant under the parity operator. 

Our interest in this paper is in the outgoing fluxes of particles N, energy E and 
angular momentum J, which means that we do not need to compute all the components 
of the renormalized stress-energy tensor. We require only the components T tr and T rip , 
and our analysis is further simplified by the fact that these components do not require 
renormalization. This latter fact is proved for a quantum scalar field in [50J. Their 
argument is, however, purely geometric, and can readily be extended to the case of an 
electromagnetic field, thereby showing that the counterterms in Ref. |Mj vanish for these 
components of the stress-energy tensor. In addition, a similar argument can be used to 
show that the linearly divergent term [52] also vanishes for these components. We may 
then use the results jSOJ EH] 



dE 
~db 
dJ 
~dt 



So 



So 



d9 dpr 2 sin 9 (U-\T tr \U-), 

d6 dip r 2 sinO (U-\TZ\U~), (27) 



where the integration is taken over the sphere at infinity. The contribution to T tr for an 
"up" mode takes the form: 

T imw f w q2 Q 2 \ , fl3 c os sin 2 9 . 2 2 , 

+asin9 (-iS/ mu 99_iS{ miJ — +iSi muJ de + iS£ muJ ) j. (28) 

In the above equation, Ti muj is the transmission coefficient, defined as a ratio of energy 
fluxes: 

with E^ nc ^ and _E( tra ) denoting, respectively, the incident energy on the black hole and 
the energy transmitted down the event horizon. The quantity T^ ma; is the analogue of the 
"absorption probability" defined in [21] . The transmission coefficient Ti mil) can be written 
in terms of -iR*) 1 ^^, the Wronskian of the "in" mode radial functions, as follows 

m-. 

— i 

T£ mw = ^4— gWt+l-R, -lR*]fmwi (3°) 

where we have chosen the normalization -i-R^™ c = -i-R"m™ C = 1- The contribution to 
T tr for an "in" mode is the same as for an "up" mode but with a change of sign. There 
are also similar expressions describing the contributions to T rip for both the "in" and "up" 
modes. 

For the fluxes in which we are interested, only the behaviour of the stress tensor at 
infinity is needed. As r — > 00, the expression (|28|) simplifies to: 



(Ttrfy^ + 'PTtr^J ~ -4^272"^ ("l^W + +lSirm>) ■ ( 31 ) 
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It is important to note that both the positive and negative helicity angular functions 
contribute to this expression. In j2J, only one helicity was included for each mode, 
which resulted in the angular distributions in that paper being highly asymmetric about 
the equator. By using the correct expression given above, we will find later that the 
angular distributions are in fact symmetric about the equator, as would be expected from 
a physical point of view. 

We now use (j31|) and the corresponding result for T rtp to find the following formulae 
for the differential emission rates of particles, energy and angular momentum, integrated 
over all angles 9, as a function of the mode energy uo: 



d 2 N 



1 x 1 

^EE^T7T T -' ( 3 2) 



dtduo 2% 

1=1 m=-£P=±l 



d 2 E 



dtduo 2n 

i=\ m =-lP=±\ 



d 2 J 



1 


exp 




- 1 




CO 




exp 


(w/Th) 


- 1 




m 




exp 




- 1 



-y X I 

— r,r,/T^_1 T ^' (33) 



dtduo 2w 

l=\ m=-£P=±l 

as well as the angular distributions of particles and energy: 
d 3 N 



d(cos6)dtduo ~ ^^^ t ^ i ex P (u/T n )-l T ^^ 1 ^ ++1 ^ 



, (35) 



d(cos 6)dtduo 4vr ^ ^ ^ exp (w/T H ) - 1 

where Tfi is the Hawking temperature of the (4 + n)-dimensional, rotating black hole, 
given by 

_ ( n +l) + (n-lK ^ 
TH " Ml + a'K • (37) 

An important part of the above expressions is the sum over the polarization P, which is 
in agreement with |2S1 • The transmission coefficient Te muJ is independent of P, so in effect 
we simply multiply by a factor of 2. As in [2H], we do not study the angular behaviour 
of T rip (corresponding to an angular distribution of the flux of angular momentum) as it 
will not be directly observable. However, for the interested reader, T rtp is studied for the 
4-dimensional Kerr black hole in |3()j . 

Let us finally note here that the emission of Hawking radiation will be assumed to 
take place through the emission of gauge particles with energy much smaller than the 
mass of the black hole. This will guarantee that there is no significant back reaction to 
the gravitational background due to the change in the black hole mass after the emission 
of a particle. Since the black hole radiation spectrum is centered around the temperature 
of the black hole, we must demand that Th Mbh- This eventually translates to 
Mbh 3> M*, a condition that has already been imposed to guarantee the absence of any 
quantum corrections to our analysis. 



4 Numerical Analysis 



The greybody factor, or transmission coefficient T^ mw , given by Eq. ([3U|h can be calculated 
from the solutions of the radial Teukolsky equation (|17|) for propagation on the brane 
black hole background; however, we still need to solve the angular equation (fTSj) since 
its eigenvalue appears in the radial equation. We will also require the spin-weighted 
spheroidal harmonics themselves for the angular distributions of the particle and energy 
fluxes (I35j l36j) . In this section, we briefly outline the numerical methods used to solve 
both the angular and radial equations of a spin-1 field in the aforementioned black hole 
background. 

In order to solve the angular equation, we use the shooting method as described in 
|33] . This method is based on the one applied in j2Hl E3j to solve the spin-0 spheroidal 
harmonic equation, and here we adapt it to the non-zero spin case. We start by using the 
Frobenius method, and re-write the angular function as follows: 



hStmvix) = (l-x) a+ (l + x) a - h y emu) (x), (38) 



where 



\m + h\ \m — h\ 



a + = , a. 



(39) 



2 ' 2 
The function hUimuj{.x) then satisfies the differential equation 

d 2 d 
x 2 )— - - 2 [a + - a_ + (a + + a_ + l)x] — + h\muj + Hh + 1) - a 2 u 2 
dx z ' dx 

+ 2mau — (ce + + a_)(a + + a_ + 1) + a 2 u 2 x 2 — 2haux\hyinw{ x ) — • 

(40) 

From Eq. pOj). it is clear that hVirw cannot be calculated at exactly the end-points x — ±1 
since these are (regular) singular points. We will, therefore, integrate this equation from 
the point x\ = — 1 + e up to the point X2 = +1 — e, where e<l. 

The shooting method for solving Eq. (jlUj) requires an ansatz, hXimu, for the eigen- 
value h^imuj of the angular equation. Imposing regularity of the solution at x = — 1 
determines the values of hVimu and hy\ m w a ^ x ii U P ^° an arbitrary overall normalization. 
The equation is then integrated from x\ to X2, as an initial value problem. If the solution 
obtained at x<i does not comply with the boundary condition of regularity at x = +1, 
the value of hXimu is improved. The differential equation is then integrated again with 
the new value of hXzmw until regularity at x — +1 is met within the desired accuracy; the 
derived solution for the spin-weighted spheroidal harmonics is then normalized so that 
Eq. (fTfij) is satisfied. 

For an ansatz h\imw for the angular eigenvalue, the numerical solution hVe^ of Eq. 
(|40p is a combination of both a regular and an irregular solution, i.e. 

hVe^Z — Mh^tmu) hVlmu + B( h \ £nuu ) hVlrrL ' (^) 

where hVtmw an d hVe^ are the regular and irregular solutions at x = ±1, respectively, 
and A and B unknown functions of hXemu- The value of hXimu must be modified so that 
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Figure 1: (a) The separation constant ^^^(aw), for the modes t = m = 1,2 and 
i = 2,m = 1, and helicity /i = — 1; (b) Spin-1 weighted spheroidal harmonics for a gauge 
field with h = —1, £ = 1, m = 0, and for aw = (0, 1, 2, 3). 



only the regular term A h yi muJ is retained. For general spin 3 , this is achieved by requiring 
that hXimuj is a zero of the function 

g{ h \^) = hv'jSZM - h y';Z IOX (x2) , (42) 

where hVl^^i^) is an approximation to the (unknown) regular value: 

«y'£TM - faL,M • (43) 

The fraction hVemui^) I , hVimu>{ x 2) can be found analytically by calculating a series expan- 
sion of hVtmw around x = +1. This series and further details can be found in [12]. As an il- 
lustrative result, in Fig. [TJ we depict the value of the separation constant hAi muJ = h^emuj + 
2mauo — a 2 u 2 , and of the spin-1 weighted spheroidal harmonics for various field modes. 

Once the eigenvalues of the angular equation ([T3*|) are found, we can proceed to solve 
the radial equation (|17|) . The radial Teukolsky equation in four dimensions contains a 
long-range potential in the spin-1 case, which cannot readily be integrated numerically. 
We use the following definition of long/short-rangeness |57j: a second-order differential 
equation 

d 2 Y{x) „, , dY(x) „. w „ , , t , 

~ix^ + A{x) + B[x) Y{x) = (44) 

is said to be short-range if, and only if, A(x) = 0(x~ d ) and B(x) = b 2 ± + 0(x~ d ), where 
d > 2 and b± are constants, when x — > ±oo. If this condition is guaranteed, then the 
asymptotic form of the solution is 

Y(x) ~ < , ., , ' (45) 



3 In the scalar case, regularity at x = +1 may be imposed by requiring that oA*mu> is a zero of the 
function g(o\£ muJ ) = oy£^{x%) ~ oyi mw (%2); this function clearly tends to zero as o-W w approaches 
the correct eigenvalue while it tends to infinity when oAf mu is far from it due to the behaviour of the 
irregular solution. This method (see also a variation in |28p can be used since, for scalar fields, the 
analytic value hy[ muJ ( x 2) is known because oy'e muJ (x) oc oy'^ muj {— x). However, for a non-zero spin field, we 
have hy'e muJ {x) oc ~hy'f mLU {~x), and therefore we do not know the analytic value h,y'i muJ { x i) f° r a particular 
value ft ^ of the helicity; therefore, an alternative method must be used. 
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Detweiler showed |55| lo^] that a particular transformation from the radial function -\Ri m u 
to a new radial function X^ muJ leads to a differential equation which is real and contains 
a short-range potential in terms of the tortoise co-ordinate r* (|24|) . This new differential 
equation, being short-range (in the homogeneous case), is suitable for numerical integra- 
tion. This method is followed and detailed in [20] for the four dimensional, Kerr-Newman 
case. 

In this work, we have followed Detweiler's method and generalized his differential 
equation valid for the Kerr black hole to any 4-dimensional line-element of the form 
including the case of a higher-dimensional black hole background induced on the brane. 
It can be easily checked that the Teukolsky-Starobinskh identities [33J EHj are valid for a 
general metric function A(r); in particular, 

+ xRim. - 2 ^- - — J ^- - — J -tRtw . (46) 

It can also be easily checked that the same is true for all the formulae in until we 
reach the point where the potential is written out explicitly, i.e. Eq. (61) in For a 
general function A, we define a new radial function X^ muJ according to 

Xe^^^ + a^A-^x^, (47) 

where 

Ximu = a -iRimu + (3 ~ 1 d ^ mUJ ■ ( 48 ) 

The functions a and (3 appearing in the definition of xtnun are quite complex. To define 
them, first write a and (3 in terms of two new quantities p and q as follows: 

a = p + a D q , /3 = b D q , (49) 

where 

a D = ~ [2K 2 + A(iK' - -iA^)] , 

UK . . 

bo = (50) 

The p and q functions are then themselves defined as 

(2K 2 ^ ~ 1/2 

P = lBimuj ( -l^imw + iBimu) 

pA 

q = V7^> (51) 

iB? muJ = -iAL, + Amau - Aa 2 u 2 . (52) 
The function Xi muj now satisfies the differential equation 



where, finally, 



d 2 X, 



dr 



— -iUXf muJ — , (53) 
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with 



-iU 



A_iA 



tmu) 



K 2 



[r 2 + a 



2\2 



+ 



A 2 (KpJ 


( A ^ 


j 3/2 


" a 1 / 2 " 


// 


(r 2 + a 2 ) 2 Kp 


\r 2 + a 2 / 




_(r 2 + a 2 )V2_ 





(54) 



Note that the signs of ap and bp above differ from those in py. so that the sign in the 
Teukolsky-Starobinskh identity (f4T)J) agrees with that in 

The explicit form of the potential (|54|) given by Detweiler uses A = r 2 — 2Mr + a 2 , 
and therefore it only applies to the four dimensional, Kerr background. Here, we have 
calculated it for a general function A(r), and the result is: 

[uj(r 2 + a 2 ) - am} 2 A_iA^ maJ 

— 1V» — / r, ^T^ r 



+ 



(r 2 + a 2 ) 2 (r 2 + a 2 ) 2 

A {2(2r 2 - a 2 ) A + (r 2 + a 2 ) [A"(r 2 + a 2 ) - 2rA']} 
2(r 2 + a 2 ) 4 

A {2(5r 2 + z/ 2 )A - (r 2 + z/ 2 ) [5rA ; + + z/ 2 )] } 

(r 2 + a 2 ) 2 [(r 2 + z/ 2 ) 2 + 77A] 
6r(r 2 + z/ 2 ) 2 A [2rA - (r 2 + z/ 2 )A'] r]A (A') 2 [2(r 2 + z/ 2 ) 2 - r]A] 



(55) 



(r 2 + a 2 ) 2 [{r 2 + v 2 ) 2 + 77 A]' 



4(r 2 + a 2 ) 2 [(r 2 + z/ 2 ) 2 + r^A] 



2 • 



where 



v 



V = 



a 2 — arajuj , 

2^ 



(56) 



When A(r) is given by Eq. (J2J), the potential _iW is found to be real and short-range, for 
any number of extra dimensions n > 0. That is, the differential equation is short- 
range, in terms of r*, in both limits r — > +00 and r — > for any n > 0. Furthermore, it 
can be checked that this differential equation is short-range, in terms of the independent 
variable r, in the limit r — > +00, for n > 1. Note, however, that it is not short-range in 
terms of r in the limit r — > for any n > 0. 

The relation r* = r*(r) can be easily found in closed form in the cases n = — ► 3 but 
not in the cases n = 4 — > 7. We, therefore, found it useful to define a new variable ta such 
that behaves like r*, for r — > r^, and like r, for r — > +00 (to ensure short-rangeness) , 
and such that we can find the relation = r&(r) in closed form Vn > 1. The differential 
equation will be short-range (for n > 1) with respect to this new radial variable, both for 
r — > and r — ->• +00. A radial variable ta that satisfies these properties can be defined 
via 

dr& r 



The relation ta 



dr r — rh 
r&(r) in closed form is then given by 



rA = r + r h log 



r - r h 



(57) 



(5? 



This relation can be inverted with the use of Lambert's W function jJH] to yield 

1 



W { — e {rA - rh)/rh ] + 1 



(59) 
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Note that ~ r*, as r — > r h , +00. If Eq. (jo^j) is re-expressed in terms of r^, the resulting 
differential equation is 

d 2 X lmu) ( r -r h )[rF^(r 2 + a 2 ) + F A (a 2 -r 2 )}dX emuJ fr 2 + a 2 \ 2 

-iU{r)A.£ muJ — U , 
(60) 



dr\ r 2 F&(r 2 + a 2 ) dr& V rF/\ 



where 



(61) 



r -r h 

is finite at the horizon. In this work, we numerically solved the differential equation (J53j) 
in the cases n = — > 3, and the differential equation ()6(jjl in the cases n = 4 — > 7. 

As in j2Hj, we solved the radial equation (either (J53J) or ffHU|) as applicable), to find 
the "in" modes, for which the function Xi muJ has the asymptotic forms 

XZ.~{ . (62) 
[e-^+AZ^e^ (r^+00). 

Our ultimate aim is to calculate the greybody factor Ti muj . To find an expression for this 
in terms of the function Xi muJ , we first note that the inverse of relation (j^Hj) is given in 
and it is equally valid for a general function A(r): 



Imu) 



n ( ( 1 + n - ; f a' + 



-1 



/" 1 oi\ a ^Xi 

[a + P) Xtmw - P 



dr 



(63) 



Using this relation, we can obtain the coefficients of the radial Teukolsky function 
at r — > +00, which correspond to the normalization chosen in Eq. (J52J) : 

? in : ref /I, ,2 J 

-l R tC = oa72L .1 • ( 64 ) 



R 111 ' 4r ; 2 

-l^mw _ ^ D in,inc 



pin,inc,i n o ' 1 «mu r>3/2| 

These expressions, together with Eq. (J3U|) . allow us to obtain the following simple expres- 
sion for the greybody factor: 



Ta™ = -|B£J 2 . (65) 
00 

By using the above result, we may then proceed to compute the fluxes of particles, energy 
and angular momentum, given in Eqs. (j32fj3Tlj) . by a brane- induced black hole. 



5 Numerical Results 

We now turn to the presentation of the exact numerical results, first, for the transmission 
coefficient and, subsequently, for the flux, power and angular momentum spectra for the 
emission of gauge fields on the brane from a higher-dimensional, rotating black hole. As 
in the case of the emission of scalar fields [2E], we wish to study the dependence of the 
various emission rates on the dimensionality of spacetime n and the angular momentum 
parameter a*. Also, the angular distribution of the emitted number of particles and energy 
will be derived - as we will see, the particular shape of the angular distribution comprises 
a clear signature not only of the emission during the spin-down phase in the life of the 
black hole but also of the type of emitted particles. Finally, the total emissivities of flux, 
power and angular momentum, as well as the superradiance effect, will be discussed, and 
more particularly their dependence on the two fundamental parameters n and a*. 
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Figure 2: Transmission coefficients for spin-1 emission on the brane from a 10-dimensional 
black hole with variable a*, for the modes (a) £ — 1, m — 0, (b) £ — l,m — — 1 and (c) 
e = l,m= 1. 



5.1 Transmission Coefficient 

We start the presentation of our results with the ones for the transmission coefficient 
IVmo;, as these follow from the numerical solution of the radial equations and (jHOJ) by 
using as input the exact eigenvalue -lXtmw- As we will see, the transmission coefficient 
bears a strong dependence on both the angular momentum parameter a* of the black hole 
and the number of additional spacelike dimensions n. In addition, its behaviour depends 
on the particular mode studied as well as on the energy regime that is investigated. 

Figure El depicts the dependence of Ti muJ on the angular momentum parameter a* 
- for simplicity, the number of extra dimensions has been fixed to n = 6. The three 
panels show the modes {£ — 1, m — 0), (£ = 1, m = —1) and (£ = 1, m — 1), and reveal 
the non-monotonic behaviour of Tg muJ as a* varies. For the {£ = l,m = 0) mode, the 
transmission coefficient first increases and then decreases, as a* increases. This behaviour 
is justified by, and reflected in, the behaviour of the gravitational potential ()55|) seen by 
the propagating gauge field: as a* starts increasing, the barrier lowers, thus giving a boost 
to Ti muJ ; as a* increases further, however, the one-peak potential changes to a double-peak 
form with the height of the peak located closer to the horizon increasing with a*, which 
leads to a suppression of the transmission coefficient. 

For the (I — l,m — —1) mode, the behaviour of Ti muJ is found to be w-dependent, 
too. For low values of lu, as a* starts increasing, the one-peak potential tends to become 
significantly more localized, leading to the observed enhancement of the transmission 
coefficient. For large values of u, the gravitational barrier changes to a well whose depth 
increases with a*, thus suppressing the value of Tf mul . The dependence on to is even 
more striking for the third depicted mode {£ — l,m— 1). This mode is a 'superradiant' 
mode, one that is characterized by a negative value of the transmission coefficient over 
the energy regime Cj = u — mQ < 0, according to Eq. ()65|) . and thus by a value of the 
reflection coefficient larger than unity (see also section 5.6). For fixed dimensionality, Fig. 
H[c) reveals that Ti muJ is greatly suppressed with a* over the superradiance regime but 
enhanced everywhere else. The behaviour of the potential again supports the observed 
behaviour: in the low-energy regime, an increase in a* either raises the barrier or deepens 

IF) 




the well that characterizes the gravitational potential, thus suppressing the transmission 
coefficient to even smaller values; in the high-energy regime, as a* increases, the potential 
is characterized by the simultaneous presence of a barrier and a well, whose height and 
depth, respectively, decrease as a* increases, thus justifying the observed enhancement. 

The dependence of the transmission coefficient on the dimensionality of spacetime 
is depicted in Fig. El for a fixed value of the angular momentum parameter, i.e. a* = 1. 
The left panel depicts the behaviour of a non-superradiant mode - in this case the £ = 1, 
m = — 1 mode - and reveals the suppression of Ti muJ , as n increases. The above monotonic 
behaviour holds for all non-superradiant modes, with m < 0, and in all energy regimes. 
The right panel of Fig. |3] shows the behaviour of the corresponding superradiant mode 
(I = 1, m — 1), in which we observe that the supression of T^ ma; with the number of extra 
spacelike dimensions holds also for this type of modes. The suppression is present also in 
the superradiant regime, where an increase in the value of n leads to a further suppres- 
sion (larger negative value) of the transmission coefficient. In terms of the gravitational 
potential, one may easily see that, for all modes, the height of the barrier - or, the depth 
of the well, depending on the values of a* and to - increases as n increases, thus justifying 
the behaviour depicted in Fig. 

The various emission spectra - particle, energy and angular momentum - as well 
as their dependence on the parameters a* and n, will be determined by the behaviour of 
both the transmission coefficient Ti mw and the thermal factor (e W//TH — appearing 
in Eqs. (J32j) -(|M |) . Having discussed the behaviour of the transmission coefficient, let us 
now comment on the dependence of the thermal factor on n and a*. The behaviour of 
the latter on the number of additional dimensions is encoded in the expression of the 
temperature of the black hole, given in Eq. (J37|) . A simple study reveals that, for a* 
and r/j fixed, dT^/dn > 0, and the temperature of the black hole always increases with 
the value of n, thus giving a boost to the thermal factor and, therefore, to all emission 
spectra. As the thermal factor behaves in exactly the opposite way from the transmission 
coefficient in terms of n, it remains to be seen which of the two is the prevailing factor 
that determines the behaviour of the final emission spectra. 

The dependence of the thermal factor on a* is more elaborate: in this case we find 



Ifi 



that, for n and r h fixed, 



d 



lo — mVt 



) 



4irr h 



mal (n — 1) + 4cl*uj — m (n + 1) 



(66) 



da* 



[( n+ l) + ( n _l) a 2]2 



according to which the behaviour of the thermal factor, as a* varies, depends on the 
particular mode, the energy and angular-momentum regime, and the dimensionality of 
spacetime. For modes with m = 0, the behaviour of the (u — mfl) /Tfi is monotonic, and 
increases with a* for all values of n and cu, leading to the suppression of the thermal factor 
and thus of the various emission spectra. Combining this result with the behaviour of the 
transmission coefficient depicted in Fig. Ufa), we may conclude that the m = modes 
will significantly contribute to the various spectra only for black holes with an angular 
momentum in the low a*-regime, where the suppression of the thermal factor is small and 
the transmission coefficient is enhanced. Turning to the modes with m < 0, we find that, 
for n < 1, the r.h.s of Eq. (jfifij) is always positive, and the various spectra are suppressed 
for all a* and u; for n > 1, the thermal factor is found to be again suppressed in the low 
a*-regime but significantly enhanced in the high a*-regime and low tu-regime. According 
to Fig. 12(b), the transmission coefficient will also be enhanced for these values of a* and 
lu, therefore, modes with m < will be mainly emitted in the low-energy regime from 
black holes with a large angular momentum. Finally, for modes with m > 0, we need to 
distinguish between the superradiant and non-sup err adiant energy regimes. In the latter 
one, we find that the thermal factor is actually suppressed with a*, while, according to 
Fig. 12(c), the transmission coefficient is enhanced; it is, therefore, the 'battle' between 
these two quantities that will define the final spectrum for modes with m > in the 
non-sup err adiant regime. For modes with m > and energy in the superradiant regime 
uj — mQ < 0, the situation changes radically: the thermal factor is found to be suppressed 
with a*, however, its value is now negative. Since this is the energy regime over which the 
transmission coefficient also acquires a negative value, the Hawking radiation flux remains 
positive and it is actually significantly enhanced with a* as both the thermal factor and 
Ti muJ are driven towards larger negative values, as a* increases. 

In the next three subsections, we present the various emission spectra for spin-1 
fields by using the results derived above. Their final expressions (|32|) - (j36|) involve a sum 
over both m G [—£,£] and I £ [1, oo), thus rendering them mode-independent. Although 
this makes the contribution of individual modes impossible to discern, the derived spectra 
will be the combined effect of all radiated modes and contributing factors. 

5.2 Flux emission spectra 

The number of gauge bosons emitted by the black hole on the brane per unit time and 
unit frequency was computed first, by using Eq. (J52J). Our results are presented in Figs. |U 
and El in terms of the dimensionless 4 energy parameter ur^, and for fixed dimensionality 
and angular momentum parameter of the black hole, respectively. 

Figures E(a,b) depict the particle emission rate for a 5-dimensional (n = 1) and 
a 10-dimensional (n = 6) black hole, respectively, and for various values of a*. Both 
plots reveal the clear enhancement in the number of particles emitted per unit time and 

4 Throughout our analysis, we will be assuming, for convenience, that the black hole horizon radius 
remains fixed, as the parameters n and a vary. 
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Figure 4: Flux emission spectra for spin-1 particles on the brane from a rotating black 
hole, for (a) n = 1, and (b) n = 6, and various values of a*. 




Figure 5: Flux emission spectra for spin-1 particles on the brane from a rotating black 
hole, with a* = 1, for various values of n. 



frequency by the black hole, as the angular momentum of the black hole increases. In both 
cases, this enhancement is observed over the whole energy band, but the exact form of the 
spectrum depends strongly on the dimensionality of spacetime: for low values of n, there 
is a clear preference for the emission of low-energy gauge bosons, although the emission of 
particles with higher energy becomes increasingly more likely, as the angular momentum 
increases; for high values of n, the emission of low-energy quanta is still favoured but the 
probability for the emission of high-energy ones has become substantially more important 
than the one in the low-dimensional case. The spectra in both cases are characterized by 
oscillations caused by the higher partial waves coming gradually into dominance as the 
energy increases. 

As is well known from previous studies, an increase in the dimensionality of space- 
time causes a significant enhancement in the number of particles emitted per unit time and 
frequency by the black hole; this effect was observed in the case of a higher- dimensional 
Schwarzschild [TTJ [JjJJ Ej , Schwarzschild-de-Sitter [T3], and Schwarzschild-Gauss-Bonnet 
[To] black hole, as well as in the case of scalar emission from a higher-dimensional Kerr-like 
black hole I2H] . The same enhancement is anticipated in the present case of the emis- 



18 




Figure 6: Angular distribution of the flux spectra for gauge boson emission on the brane 
from a rotating black hole, for n — 1 and a* = (0, 0.6, 1). 




Figure 7: Angular distribution of the flux spectra for gauge boson emission on the brane 
from a rotating black hole, for n = 6 and a* = (0, 0.6, 1). 



sion of gauge bosons, and this behaviour can be readily seen by comparing the vertical 
axes of Figs. Ufa) and (b). A more complete picture can be drawn from Fig. where 
the flux spectrum is given for fixed angular momentum (a* = 1) and variable n. As we 
see, the particle emission curves become significantly taller and broader, as n increases, 
which signals the enhancement of the total emissivity of the black hole by virtually orders 
of magnitude. 

We now turn to the angular distribution of the flux spectra: we study again the 
two indicative cases of n = 1 and n = 6, in an attempt to detect potential changes in 
the angular distribution of particles as the dimensionality of spacetime changes. In each 
case, we present spectra for three different values of the angular momentum parameter, 
i.e. a* = (0, 0.6, 1); this will help us see more clearly the effect of the black hole rotation 
on the angular distribution pattern. These spectra are shown in Figs. El and [H For both 
values of n, the spectra show no angular variation for a* = 0, as expected. However, as 
a* starts increasing, a much more interesting pattern emerges: low-energy gauge particles 
tend to be emitted along the rotation axis {6 = ±7r), but this tendency soon dies out as 
the energy of the particles increases. The concentration of the emitted gauge particles 
with low energy close to the rotation axis is caused by the non-trivial coupling between 
the spin of the particles and the angular momentum of the black hole, a coupling that 
was absent in the case of emission of scalar particles [2H] ■ For higher-energy particles, this 
effect gradually becomes sub-dominant compared to the centrifugal force generated by the 
rotating black hole; this forces the emitted particles to concentrate around the equatorial 
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Figure 8: Power emission spectra for spin-1 particles on the brane from a rotating black 
hole, for (a) n = 1, and (b) n = 6, and various values of a*. 




region (8 = 7r/2) instead, an effect that seems to be enhanced as both the dimensionality 
of spacetime and angular momentum of the black hole increase. 



5.3 Power emission spectra 

We now proceed to the power spectrum, i.e. the energy emitted by the black hole in the 
form of gauge bosons per unit time and unit frequency. This is given by Eq. (JHHj) . and it 
is displayed in Figs. [H] and El for fixed dimensionality and angular momentum parameter, 
respectively, as in the previous subsection. 

Figures|Hf a) and (b) depict the energy emission spectra for a 5-dimensional and a 10- 
dimensional black hole, respectively, and for variable a*. In both cases, the enhancement of 
the emission rate, with the value of the angular momentum parameter, is obvious over all 
energy regimes 5 . In the n = 1 case, the localized - around the low-energy regime- curve, 

5 The same enhancement in the emission rate of gauge particles was found in the low-energy regime in 
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Figure 10: Angular distribution of the power spectra for gauge boson emission from 
rotating black holes, for n — 1 and a* = (0, 0.6, 1). 




Figure 11: Angular distribution of the power spectra for gauge boson emission from 
rotating black holes, for n = 6 and a* = (0, 0.6, 1). 



for a* = 0, gives its place to a broad, oscillatory curve, for higher values of a*. Although 
the emission rate is gradually decreasing, as the frequency of the particle increases, it 
retains more than 50% of its value when the frequency has increased by a factor of 5. In 
the n = 6 case, the situation has changed even more radically: a peak around the low- 
energy regime can still be seen in the spectrum - remnant of the dominant peak appearing 
in the corresponding flux spectrum; however, it is the emission of gauge particles with 
much higher energy that is the most effective emission channel now. 

The additional enhancement, caused by the increase in the dimensionality of space- 
time, can be easily seen again by comparing the axes of Figs. Ufa) and (b). This enhance- 
ment is shown in detail in Fig. |§1 for fixed a* and n = 1, 7. The increase in both the 
height and width of the emission curves, as the number of additional spacelike dimensions 
increases, can be clearly seen. From the above, we anticipate a strong enhancement of 
the total emissivity, in the form of gauge particles, of a higher-dimensional rotating black 
hole with parameter n, a point to be studied in more detail in subsection 5.5. 

The angular distribution of the power spectra for gauge boson emission on the 
brane is shown in Figs. El and ^fl for the cases of n — 1 and n = 6, respectively, 
and for a* = (0,0.6,1). The same features exhibited by the flux spectra appear also 
here: the angular variation - absent for a* = - makes its appearance in the spectrum 

|21) by following an analytical approach valid in the low-energy and low- angular- momentum limit. The 
suppression shown in Figs. 6 and 7 of |21j over the intermediate and high-energy regimes are an artificial 
result only, caused by the breakdown of the approximations used, and therefore should not be trusted. 
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Figure 12: Angular momentum spectra for spin-1 particles from a rotating black hole, for 
(a) n = 1, and (b) n — 6, and various values of a*. 



as soon as the angular momentum parameter acquires non-zero values. Again, the low- 
energy emission channel runs parallel to the rotation axis, while, for higher-energies, the 
emission concentrates in the equatorial region. The amount of energy emitted per unit 
time and frequency along the equatorial plane increases also with the angular momentum 
parameter and the dimensionality of spacetime. 

We should note here that our spectra are symmetric under the change 8 — > n — 8, 
contrary to what was found in [21]; the difference is due to the fact that, in our analysis, 
both polarisations of the gauge field have been taken into account. Although, the two 
polarisations may have a different coupling with the black hole rotation, these couplings 
are both reversed when moving from one hemisphere to the other, thus resulting in a 
symmetric spectrum [30] . 

5.4 Angular momentum spectra 

The final emission rate that we would like to study is the one of the angular momentum 
of the black hole. This differential rate, per unit time and unit frequency, is first shown in 
Figs. ITW a.b) for n = 1 and n = 6, respectively, and for variable a*. Its behaviour is very 
similar to the one exhibited by the power emission spectrum. For n = 1, the process of 
the angular momentum loss takes place predominantly through the emission of low-energy 
gauge particles, although the emission of particles with higher energy is still substantial. 
For n = 6, on the other hand, the black hole loses its angular momentum through a 
more democratic process of emission, with a slight preference towards the emission of 
high-energy particles. 

Figure IT31 depicts the dependence of the angular momentum loss rate on the dimen- 
sionality of spacetime, while fixing the angular momentum parameter at a* = 1. As in 
the previous two cases of flux and power spectra, this rate is greatly enhanced, as the 
dimensionality of spacetime increases. The exact amount of the enhancement of the total 
angular momentum emissivity - which, as we will see, can reach more than an order of 
magnitude - will be discussed in the next subsection. 



77 



° 0.5 1 1.5 2 2.5 3 3.5 4 

cor h 

Figure 13: Angular momentum spectra for spin-1 particles from a rotating black hole, for 
a* = 1 and various values of n. 



5.5 Total emissivities 

The total emissivities follow by integrating the various differential emission rates over the 
frequency regime, and give the total emission rate by the black hole per unit time. In 
principle, the integration over the frequency must cover all regimes in which the black hole 
radiates. However, deriving the whole Gaussian emission curve, as the angular momentum 
parameter and dimensionality of spacetime increase, soon becomes an unrealistic task. By 
looking, for instance, at Figs. 03 E] and El one may easily see that while, for low values of 
n, all emission curves have been completed by the time the value of the energy parameter 
ujrh = 4 is reached, for high values of n, the emission curves are still very close to their 
peak value. As both a* and n take larger values, more and more high-energy particles are 
emitted by the black hole, and the tail of the curve extends to values of urh much larger 
than the cutoff value of 4, up to which the radial part of the equation of motion of the 
gauge field has been numerically solved. Extending our energy regime towards much larger 
values has proved to be an extremely time-consuming process, therefore, in what follows, 
we will derive the "total emissivities" by integrating the emission rates up to our chosen 
cutoff value. As is obvious from the results presented in the previous subsections, an 
increase in the two fundamental parameters (a* and n) leads to a significant enhancement 
of the emission rates over the whole energy regime; the high-energy regime left out of our 
computation of the total emissivities, for large values of a* and n, would only increase 
further this enhancement. Therefore, the "total emissivities" derived, and presented in 
this subsection, can be considered as robust lower bounds to the exact ones. 

By following the above method, we have computed the total emissivities for the flux, 
power and angular momentum emitted by a higher-dimensional, rotating black hole on 
the brane in the form of gauge fields. The results, for fixed angular momentum parameter 
(a* = 1) and dimensionality (n = 6), are depicted in the two histograms in Figs. [T4T a) 
and (b), respectively. The enhancement of all three emissivities, that was foreseen while 
commenting on the results derived in the previous subsections, can be clearly observed. 
In the first histogram, the value of the number of additional spacelike dimensions n varies 
from 1 to 7. A simple calculation shows that the number of particles emitted per unit 
time by a 11-dimensional, rotating (with a* = 1) black hole is 30 times larger than for 
a 5-dimensional one. Similarly, the total amount of energy emitted per unit time and 
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Figure 14: Total emissivities for spin-1 particle emission on the brane from a rotating 
black hole as a function of (a) n, for a* = 1, and (b) a*, for n = 6. 



the angular momentum loss rate for a 11-dimensional black hole are more than 50 and 
30 times, respectively, larger than for a 5-dimensional black hole with the same angular 
velocity. For the reason explained in the beginning of this subsection, we expect these 
numbers to increase considerably when the high-energy part of the spectrum, for n = 7, 
is included in the calculation. 

The second histogram in Fig. [J3 depicts the enhancement of the three emissivities 
when the dimensionality of spacetime is kept fixed (at n = 6) and the angular momentum 
parameter varies. The data, used to produce this plot, reveal that the number of particles 
and amount of energy emitted by a 10-dimensional black hole with a* = 1 is 4 and 5 
times, respectively, larger than the ones for a 10-dimensional, non-rotating black hole. 
As, in the case of n = 6, a substantial part of the spectrum has been left out of the 
computation of the total emissivities, the aforementioned numbers should be treated as 
lower bounds only, with the exact numbers expected to reach values greater than these by 
even orders of magnitude. This can be deduced from the fact that, in the case of n — 1, 
where all parts of the spectrum have been successfully derived and thus included in the 
computation of the exact total emissivities, the number of particles and amount of energy 
emitted per unit time by the 5-dimensional black hole with a* = 1 is already 10 and 15 
times, respectively, larger than for its non-rotating, 5-dimensional analogue. 

From the first histogram in Fig. El we observe that the emission rate of angular 
momentum remains always larger than the energy emission rate, with the former becoming 
increasingly larger than the latter as n increases. This seems to support the argument 
that, during its decay, the angular momentum of a higher- dimensional black hole will be 
radiated away faster than its mass, thus signalling the existence of a subsequent non- 
rotating phase in the life of the black hole. Nevertheless, the second histogram in Fig. 
ITU shows that, for fixed dimensionality, the angular momentum emission rate strongly 
depends on the angular momentum of the black hole: while this rate is again significantly 
larger than the energy emission rate for rapidly rotating black holes, it is significantly 
smaller for slowly rotating ones. We might, thus, conclude that the evolution of rotating 
black holes goes through two different back-to-back stages: an initial one, where the black 
hole loses predominantly rotational kinetic energy, and a second one, where the - now 
slowly - rotating black hole loses mainly its mass. The exact duration of the two phases, 
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Figure 15: Power emissivity for spin-1 emission on the brane from a rotating black hole 
as a function of cosO, (a) for n = 6 and variable a*, and (b) for a* = 1 and variable n. 



as well as of the subsequent non- rotating one - if existent, will strongly depend on the 
initial angular momentum of the black hole and the dimensionality of spacetime. 

We finally present the integral over the frequency of the angular distribution of the 
power spectrum (J36|) to derive the emissivity solely as a function of the azimuthal angle 9. 
This is depicted in Figs. IToT a.b). as a function of cos 9 for fixed n and variable a*, and vice 
versa. As noted before, the shape of the angular distribution pattern is formed under the 
effect of the spin-rotation coupling and the centrifugal force, with the former dominating 
for low a* and low n, and the latter dominating for high a* and high n. The two symmetric 
peaks, due to the spin-rotation coupling, around cos 9 = ±1 cannot be seen in Figs. fTBIdue 
to their very small magnitude - they barely survive the integration over the frequency due 
to the very small energies to which they correspond. Nevertheless, the aforementioned 
balance between the two competing forces leads to the type of distributions shown in 
Figs. EE where a twin-peak pattern is formed. The derived shape of the integral over the 
frequency of the angular distribution of the power spectrum agrees with the corresponding 
one derived in the case of a 4-dimensional, rotating black hole [3Tij . 

5.6 Superradiance 

In this final subsection, we address the issue of superradiance j^H], i.e. the amplification 
of the amplitude of an incident bosonic wave by a rotating black hole. This takes place 
when the transmission coefficient becomes negative over a particular energy regime, i.e. 
when uo — mQ < 0, thus allowing for values of the reflection coefficient, defined as M.e muJ = 
1 — T imuJ , larger than unity. As noted in section 5.1, this behaviour coincides with the 
denominator in the expressions for the various fluxes (|32j) - (J34j) becoming also negative, 
therefore the Hawking radiation emission rates for these modes remain positive. The 
superradiance effect is known to take place in the pure 4-dimensional case for scalars, 
gauge bosons and gravitons [3B1 EE EH! ■ I n the case of a higher-dimensional black hole, 
the superradiance effect associated with either bulk [20] or brane [221 123 121] scalar fields 
has also been investigated. In [22], this effect was studied in detail in terms of both 
the dimensionality of spacetime and the angular momentum of the black hole. It was 
found that the energy amplification of 0.3% for scalars, for a 4-dimensional maximally 

9.5 




Figure 16: (a) Superradiant scattering of gauge bosons by a maximally rotating (a = 
0.99999Mb#) 4- dimensional black hole; (b) Super-radiant scattering of gauge bosons on 
the brane by an 11-dimensional black hole with a* = 1. 



rotating black hole with a* ~ 1 [60; , increased by more than an order of magnitude for a 
6-dimensional black hole, with the same angular momentum, emitting scalar fields on the 
brane. The amplification reached its highest value for the maximum values of both n and 
a* considered in that work, namely n = 6 and a* = 4, and it was found to be around 9%. 

By following the conventions of Ref. [SHI, m Figs- fTBT a) and (b), we depict the 
energy amplification of spin-1 fields due to the superradiance effect, for the cases of n = 
and a* = 0.9955, and n = 7 and a* = 1, respectively 6 . The vertical axis in these plots 
is minus the transmission coefficient, — Ti muJ1 expressed as a percentage. The horizontal 
axis is the frequency interval over which T^ ma; is negative, divided by mfl to make the 
superradiant regime (0,1) for all modes. Figure IToT a) shows excellent agreement with the 
results produced for n = and a = 0.99999M## in [3S] : the maximum amplification 
occurs for i — m — 1, at uj/VL = 0.88, and it has a magnitude of 4.4%. However, as 
the dimensionality of spacetime increases, even while the angular momentum remains the 
same, these numbers change considerably: according to the data plotted in Fig. ITBT b). 
the maximum amplification occurs again for the mode I — m — 1 but this is now at 
uj/Vt ~ 0.8, and it reaches the figure of 21.7%. The behaviour found in section 5.1 - and 
depicted in Fig. Efb) - can easily explain the enhancement with n in the superradiance 
amplification: for a* fixed, the transmission coefficient Te muJ was found to take increasingly 
larger negative values in the superradiance energy regime, as n increased, due to the 
increase in the height of the gravitational barrier. 

As the superradiance effect is clearly associated with the rotation velocity of the 
black hole, we expect its magnitude to increase with the angular momentum parameter. 
Our findings in section 5.1 support this expectation as they clearly revealed the supression, 

The value a* = 1 (equivalent to a = Mbh) corresponds to a maximally rotating black hole in the 
4-dimensional case. This extremal case is highly unstable from the point of view of numerical analysis, 
therefore a close enough value, a = 0. 99999Mb//, equivalent to a* — 0.9955, has been used in our analysis, 
as was also done in |38j . 
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Figure 17: Transmission coefficients for spin-1 emission on the brane, for n — 6, from a 
black hole with (a) a* = 0, (b) a* = 0.9, and (c) a* = 1.2. The curves, from left to right, 
correspond to the modes I = m = 1, 2, 3, 4, and 5. 



with a*, of the transmission coefficient, leading to increasingly larger negative values (see 
Fig. El). Figure IT7I depicts the behaviour of the transmission coefficient Ti mu} for a whole 
set of spin-1 modes, for a 10- dimensional black hole (n = 6) and various values of a*. For 
the case of a* = 0, shown in Fig. IT7lf a). the transmission coefficient takes only positive 
values, and no superradiance occurs as expected. For a* = 0.9, shown in Fig. fTTT b) . the 
transmission coefficient does acquire negative values over a particular energy regime before 
it crosses the zero value and passes into the positive values regime: in this case, as it can 
be seen from the different curves, the maximum amplification takes place for I = m = 1, 
and it amounts to 15.2%. As the angular momentum parameter increases further, the 
largest negative value of the transmission coefficient increases, too: for example, for the 
case of a* = 1.2, shown in Fig. ITTT c). the superradiance amplification reaches the figure 
of 36.5%. 

As a* and/or n increase further, we expect the superradiance effect to become even 
more important. In an attempt to determine the maximum value of the energy amplifica- 
tion due to the superradiance, the extreme case of n = 7 and a* = a™ ax = 4.5, as dictated 
by Eq. (J7j), was studied. Contrary to the case of the energy amplification of scalar fields 
|2*2*] . it was found that even for large values of n and a*, the lowest superradiant mode 
t = m = 1 is still the dominant one. For this particular mode, it was found that the 
energy amplification reaches the amazing figure of 4200%, at oj/Vt = 0.84. 



6 Conclusions 

The potential detection of Hawking radiation from a higher- dimensional black hole formed 
in the context of a brane-world model (either with Large Extra Dimensions PQ, or warped 
extra dimensions jHIj in the limit of d-Ads, with iAds the AdS radius) demands the 
derivation of exact results for the various emission rates during both the spin-down and 
Schwarzschild phases in the life of the black hole. After the Schwarzschild phase was 
exhaustively studied, the spin-down phase eventually came into the foreground, with a 
few early and recent papers offering partial studies of the emission of Hawking radiation 
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from a rotating brane-world black hole. The first complete study of the emission of scalar 
fields on the brane by such a black hole was performed only recently in j2H] • The present 
paper complements and expands the latter one by focusing on the emission of gauge 
bosons from a (4 + n)-dimensional rotating black hole on the brane - where the gauge 
bosons are restricted to live by assumption. A comprehensive analysis of the different 
emission rates, in terms of the values of the fundamental parameters of the model, was 
offered, and additional important aspects of physical processes associated with a rotating 
black hole, such as the angular distribution of the emitted particles and the superradiance 
effect, were also studied. 

By generalizing techniques that were originally developed for the case of 4-dimen- 
sional black holes, to apply in the cases of higher- dimensional black hole line-elements 
projected onto a brane, we first derived the equation of motion for spin-1 particles prop- 
agating in the induced-on-the-brane gravitational background, and decoupled the radial 
and angular part (apart from the presence of the angular eigenvalue that appears in both 
parts thus linking them). We then proceeded to review, and generalize where necessary, 
the principles leading to the expressions of the transmission coefficient and of the var- 
ious Hawking radiation emission fluxes. The calculation of the transmission coefficient 
demanded the numerical integration of the radial part of the equation of motion, while 
the angular distribution of the emitted particles demanded the numerical integration of 
the angular part - before, however, these two tasks could be performed, the exact value of 
the angular eigenvalue had to be computed, again numerically. Our numerical techniques 
were discussed in detail in Section 4. 

In Section 5, we presented exact numerical results, first, for the transmission co- 
efficient and, subsequently, for the particle, energy and angular momentum fluxes. The 
differential emission rates per unit time and frequency, integrated over the azimuthal angle 
6, were computed first, and their dependence on the dimensionality of spacetime, the an- 
gular momentum of the black hole and the energy of the emitted particle was investigated. 
It was demonstrated that all three emission rates were enhanced over the whole energy 
regime as the number of additional spacelike dimensions and the angular momentum of 
the black hole increase. In the case of the particle flux spectra, it was found that particles 
with a low energy are more likely to be emitted than particles with intermediate and high 
energy, however the emission of the latter ones becomes increasingly more significant as 
either n or a* increase. In the case of power and angular-momentum emission spectra, for 
low values of n, both of these quantities are predominantly reduced through the emission 
of low-energy particles; as the dimensionality of spacetime takes larger values however, 
it is the intermediate and high-energy particles that carry away most of the energy and 
angular momentum of the black hole. 

The angular distribution of the emitted particles and energy was also studied in 
detail. That became possible after the exact values of the spin-1 weighted spheroidal 
harmonics, valid for arbitrary values of the angular momentum of the black hole and 
energy of the emitted particle, were determined. The angular variation in the pattern 
of the emitted radiation can be a distinctive signature of emission from a rotating black 
hole, but also of the spin and energy of the particles emitted. In Figs. flSlf a) and (b), 
we present for comparison the angular distribution of the energy emission spectrum for a 
6-dimensional (n = 2) rotating black hole emitting radiation on the brane, in the form of 
scalar and gauge particles, respectively. As expected, the uniform distribution of particles 
seen in the case of a non-rotating spherically-symmetric black hole is now replaced by a 
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Figure 18: Angular distribution of the power spectra, for (a) scalar particle emission, 
and (b) gauge particle emission, from a 6-dimensional black hole with a* = 1 on the 
brane. 



clearly oriented one, however the exact pattern is strongly dependent on both the spin 
and energy of the emitted particle. In the case of scalar particles [2E]> the centrifugal 
force exerted by the rotating black hole forces the scalar particles to concentrate on the 
equatorial plane (6 = 0); nevertheless, this angular variation disappears in the low energy 
part of the spectrum where the spherically-symmetric I = m = modes dominate, while it 
strongly appears in the intermediate and high-energy part of the spectrum where the non- 
spherically-symmetric modes become important. In the case of gauge particles, as we saw 
in this work, the angular distribution comes as the result of two forces: the centrifugal 
one, and the spin-rotation coupling due to the non-vanishing spin of the particle. It 
turns out that the latter one dominates for low-energy particles, forcing them to align 
parallel to the rotation axis (6 = ±7r), while it becomes subdominant compared to the 
centrifugal force in the high-energy regime, causing the particles to concentrate again in 
the equatorial plane. A similar behaviour for the emission of gauge particles was seen 
also in the 4-dimensional case |3U]; in the present case, it was found that the number of 
additional spacelike dimensions work towards enhancing further the centrifugal force. 

An important quantity associated to the emission of Hawking radiation from a 
black hole is the total emissivity, that describes the emission of either particles, energy or 
angular momentum per unit time, over the whole frequency band. In fact, in our analysis, 
lower bounds for the total emissivities were computed instead of the actual total ones, 
as the differential emission rates were integrated over the energy from zero to the chosen 
cut-off value of 4. This task was performed for all three emission rates, namely for the 
particle, energy and angular momentum fluxes. Their dependence on the azimuthal angle, 
the angular momentum of the black hole and the number of additional dimensions was 
investigated. In terms of the first quantity, both the particle and energy emissivities were 
found to have a twin-peak pattern under the influence of the centrifugal force and the 
spin-rotation coupling. In terms of the other two parameters, all emissivities were found 
to increase substantially as n and a* increased; for instance, the energy emission per unit 
time increased 5-fold, for a 10-dimensional black hole, as a* increased from to 1, and 
50-fold, for a black hole with a* = 1, when n increased from 1 to 7. 

An interesting question that arises when one combines results on the various emission 
rates for different types of fields (i.e. scalars, gauge bosons and fermions) is which type 
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of fields is preferably emitted by the black hole as both n and a* vary. In the case of 
non-rotating black holes (a* = 0), it was found that a 4-dimensional black hole 
emits mainly scalar particles, a 6-dimensional black hole prefers a democratic distribution 
of particles, while the spectrum of a 10-dimensional black hole is dominated by gauge 
bosons. A similar behaviour is expected in the case of a rotating black hole. Looking at 
the vertical axes of Figs. ITHT a.b). one can easily see that the emission rates for scalars and 
gauge bosons differ significantly, for the same values of n and a*. In order to investigate 
this in more detail, we have displayed in Figs. ITW a) and (b) the emission rates for scalars 
and gauge bosons, for fixed n (n = 1 and n = 4, respectively), and for similar values of a*. 
Figure EJa) confirms the fact that for a non-rotating black hole living in a spacetime with 
a low value of n, the emission rate for gauge bosons becomes comparable to that of scalars 
only when both polarizations are included 7 . When, however, the angular momentum of 
the black hole is turned on, the emission rate for gauge bosons becomes significantly larger 
than the one for scalar fields over the whole energy regime. Moving to the n = 4 case, 
we see that indeed, for a* = 0, the gauge bosons already clearly dominate over the scalar 
fields; when the angular momentum is added on top, the increase becomes even larger. We 
may therefore conclude that a rotating black hole strongly prefers to emit gauge bosons 
over scalar particles for all values of n. 

The final issue addressed in this work was the superradiance effect associated to the 
amplification of a bosonic wave incident on a rotating black hole. We have found that, 
in the case of gauge fields, similarly to the scalar case, this effect is present and greatly 
enhanced by the addition of extra spacelike dimensions transverse to the brane. The effect 
in this case is actually more important than in the scalar one: the superradiance effect 
increases from 4.4% to 21.7% as the dimensionality of spacetime goes from D = 4 to 
D = 11, even as the angular momentum of the black hole remains the same. Vice versa, 
for fixed dimensionality, i.e. n — 6, the energy amplification varied from 0%, for a* = 0, 
as expected, to 36.5%, for a* = 1.2. As the superradiance effect is clearly augmented by 
both n and a*, the maximum energy amplification was sought by looking at the extreme 
case of n = 7 and a* = a™ ax = 4.5: the lowest superradiant mode I — m — 1 provided 
again the maximum amplification that was found to reach the amazing figure of 4200%. 

7 In the results presented refer to a single polarization of the emitted gauge field. 
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The study of the emission of gauge fields on the brane by a higher-dimensional 
rotating black hole, performed in this work, has revealed a number of distinct features 
associated with the behaviour and magnitude of the different emission rates, the angu- 
lar distribution of particles and energy, the relative enhancement compared to the scalar 
fields, and the magnitude of the superradiance effect. Apart from the theoretical inter- 
est in studying modifications of the properties of black holes submerged into a higher- 
dimensional space, we expect these features to comprise clear signatures of the emission 
of Hawking radiation of a brane-world black hole during its spin-down phase upon suc- 
cessful detection of this effect during an experiment. Having starting with the study of 
the emission of salar fields on the brane by such a black hole [2E]> we have offered in this 
work a similar, comprehensive analysis of the case of gauge fields; we hope to return soon 
with complete results for the emission of fermions that will eventually complete the study 
of the decay of a small, higher-dimensional rotating black hole through the emission of 
Standard Model fields on the brane. 
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